The gap maximum of anisotropic superconductors 



R. Combescot 

Laboratoire de Physique Statistique de VEcole Normale Superieure, associe au CNRS et aux Universites Paris 6 et Paris 7, 

24 rue Lhomond, 75005 Paris , France 
(Received February 1, 2008) 

For a completely general anisotropic order parameter (including changes of sign), we show that 
weak coupling theory is incompatible with high values of the maximum Am of the zero temper- 
ature gap as compared to the critical temperature T c , such as those found experimentally in 
Bi2Sr2CaCu20s+s where 2 Am / T c ~ 7.5 . This gives evidence for strong coupling effects. 
In particular this comes as a major support for a spin fluctuation mechanism with strong coupling, 
if one assumes that only a repulsive pairing interaction is at work in high T c superconductors. 
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PACS numbers: 74.20.Fg, 74.25. Jb, 74.72.Bk 

The last few years have seen a very important progress in the identification of the order parameter in some high T c 
superconducting compounds. Indeed there are now quite firm experimental evidences that, in Y Ba2Cu^O^ (YBCO), 
£j ■ it changes sign and that there are nodes in the gapjlj . More recently a linear dependence of the penetration depth has 
been observed in Bi^SriCaCu^O^s (BSCCO) yj3|]7and phase sensitive experiments have also given a positive answer 
There is also clear indication from Raman scattering experiments || of nodes in the gap in Hg compounds. Finally 
a spontaneous half magnetic flux quantum has been observed quite recently in Tl^BaiCuOs+s giving evidence for 
^ . a change of sign of the order parameter. These features of the order parameter are a clear indication that there is 
some important repulsive contribution in the pairing interaction. The most obvious origin for such a contribution is 
Coulomb repulsion. 



\ Yet this answer does not provide a complete physical picture for this repulsive component. Indeed |7j this repulsion 
c"j . can appear in a direct way, or it can also be the microscopic origin of low energy antiferromagnetic fluctuations, with 
pairing mostly due to exchange of these fluctuations between electrons. A qualitative difference between these two 
mechanisms is the characteristic energy of the pairing interaction. If we are dealing with direct Coulomb repulsion, 
the typical energy entering the interaction is of order of Coulomb interaction itself, that is typically a few eV. Since the 
t— I , critical temperature and the gap are quite small compared to this energy, the pairing interaction can be considered as 
instantaneous. This implies that pairing can be very well described by weak coupling BCS theory. On the other hand 
m if pairing is due to the exchange of spin fluctuations, our characteristic energy is of the order of a spin fluctuation 
frequency, which is a few tenths of eV at most. In this case the critical temperature and the gap are no longer small 
compared to this energy, and the pairing interaction can not be considered as instantaneous. This means in particular 
that pairing has to be described by a strong coupling generalization of BCS theory [|| . Therefore we can obtain an 
' indication on the kind of repulsive interaction we are dealing with by checking, as well as we can, if the superconductor 
Q\ , is satisfactorily described by weak coupling theory or if there is a need for strong coupling effects. 
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As it is well known the consequences of weak coupling theory are much more restrictive than those of strong 
coupling. Therefore in this paper we will consider some consequence of weak coupling theory and see if it can be 
made to agree with experiments. Specifically we will deal with the zero temperature gap to critical temperature ratio. 
Naturally it is well known that, for isotropic pairing, this ratio is given by the famous BCS value 1.76 . However in 
this paper we will consider the much more general case of anisotropic weak coupling BCS theory for which no such a 
simple result exists. Actually, as we already mentionned, we can have nodes in the order parameter and in this case 
there is, strictly speaking, no gap at all for the whole excitation spectrum. On the other hand for a fixed value of the 
. ,_| \ wavevector k of the excitation, we have a gap | Ak| for the excitation energy where Ak is the order parameter. We are 
interested in the maximum Am of this gap over the Fermi surface, and its ratio Am / T c to the critical temperature. 
Our reason for investigating this ratio is that we have some good experimental data for it. Surprisingly this is not so 
much for YBCO, the most investigated high T c superconductor, where the data are not very clear although 2 Am / 
T c seems to range from 6 to 8 in most measurements. The clearest data are perhaps found in BSSCO where tunneling 
experiments perpendicular to the c axis |^| give a fairly sharp peak around 30 meV leading to 2 Am / T c « 7.5 . 
The sharpness of the peak makes unlikely a shift to higher energy due to broadening. Angular dependent tunneling 
experiments in the a-b plane ]io| give even a gap maximum reaching 40 meV. One may still worry that tunneling sees 
only a surface feature. However Raman scattering clearly samples the bulk, and it gives |Tl| , [l2]| a maximum around 
550 cm" 1 , which seems to confirm 2 Am / T c « 7.5 , although the data are broad which makes the interpretation less 
secure. Taken together these experiments are suggestive of a fairly high ratio. Raman data in Hg-1212 give similar 
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results ||. Facing these experimental data, it is of interest to investigate if they can be explained by taking into 
account the anisotropy of the gap within weak coupling theory. This is the purpose of the present paper. 

In weak coupling theory the order parameter A k is obtained by solving the gap equation : 

Wc /2ttT 

A k = / d 2 k'V k , k ,A k ,27rT ]T (u£ + |A k f r 1/2 (1) 

•* m=0 

Here the integration <i 2 k' over the Fermi surface is weighted by the local density of states [(27r) 3 Vk']^ 1 ■ The summation 
over the Matsubara frequencies uj m = {2m + X)itT is limited to a cut-off w c large compared to the maximum Am of 
| At | (and therefore to T c ). We make no assumption on the effective interaction Vk.k' so the situation we consider is 
completely general. We note also that a multiband model can be considered as a particular case of gap anisotropy so 
this kind of situation is included in our study. 

At the critical temperature T c , this equation becomes linear. We call Ao, k the normalized eigenvector of V k , k ' 
corresponding to the largest eigenvalue Ao- It gives the shape of the gap at T c . Making use of 2irT c £ l/|w m | = In ( 
1.13 w c / T c ) = 1/Ao, (valid in the weak coupling limit of large lo c / T c ), we have : 

A , k = ln(1.13^) J d 2 k'V k , k ,A , k , (2) 

Below T c , A k is obtained from Eq.(l) . Now an essential feature of this equation is that the sum over Matsubara 
frequencies is dominated by the terms uj m » Am ■ This is seen by rewriting it as : 

/w c /2ttT w c /27tT c 
d 2 k'V k , k ,A k , + /d 2 k'V k , k -A k ,[27rT « + |A k <|T 1/2 ~ ]T lm + 1/21- 1 ] (3) 

m=0 m=0 

In the second term of the right-hand side we can let uj c go to infinity because the result is convergent. Since in weak 
coupling ln(1.13w c /T c ) is large, we see that the first term dominates over the second one. Therefore to lowest order 
the shape of the gap below T c is still given by Ao. k . However the size of the gap is fixed by the second term which 
is non linear. We can obtain a still exact equation for this size by multiplying Eq.(3) by AJ k and integrating over k 
(which takes into account the local density of states), leading to : 

/u c /27rT w c /27rT c 
d 2 kA;(k)A(k)[7rT J2 « + |A k ,| 2 r 1/2 - J2 |2m + l|- 1 ]=0 (4) 

m=0 m=0 

where we have made use of Eq.(2) to eliminate the interaction (using the fact that it is hermitian). Since this equation 
does not contain large terms anymore, we can replace Ao, k by A k to lowest order. This equation can also be rewritten 
as : 

ln(^) f d 2 k|A k | 2 = f d 2 k|A k | 2 2^T J2 [KP 1 - K + |A k | 2 )- 1/2 ] (5) 

J J m=0 

In particular we obtain at T = : 

Am > Jd 2 k^(k)ln(J(k)) 

ln( x^y) = Jtfl (6) 

where we have introduced <5(k) = |A k | / Am which is the absolute value of the gap normalized to its maximum value. 
This equation has already been essentially obtained by Pokrovskii [p^| . Provided that we know the shape <5(k) of the 
gap, it gives us the maximum of the zero temperature gap Am compared to the standard BCS value 1.76 T c . We 
see that the result is not sensitive to the detailed structure of the gap since the logarithm is a smooth function. 

The above result has been obtained within weak coupling theory where the parameter uj c / T c is large. If we want 
to improve on this result, we have to consider that co c / T c is not large anymore which implies to go to strong coupling 
theory anyway. There is no way to improve consistently on this result within weak coupling theory. Nevertheless we 
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might worry that the above result is a poor approximation because the dominant term is only logarithmically large. 
Fortunately the general situation is much better. This can be seen by rewriting the exact Eq.(3) at T = as : 

A k = ln(1.13^) y"d 2 k'V k , k ,A k , + J d 2 k'V k ^A k , H^T^f) (?) 

and projecting it on the complete set of normalized eigenvectors A m k of Vk,k'- If the corresponding eigenvalues are 
A m , this gives : 

( J- - 1) y d 2 kA; il!k A k = J d 2 kA; >k A k In(ij^) (8) 

where we can think of evaluating the right-hand side to lowest order by replacing A k by Ao, k . From this equation, 
the components a m = J A* k A k of the gap on the eigenvectors with m / are small because the eigenvalues A m 
are small ( this is the weak coupling limit ) . In addition the subdominant order parameters A m will correspond in the 
general case to critical temperatures much smaller than T c , which implies A T „ << Ao for m ^= ( a specific example 
of this can be found in the recent work of Palumbo et al. |14] where, within a given channel, Ai sw 0.1 Ao is found 
for all the channels). This implies that a m is reduced by a factor A m / Ao with respect to a naive evaluation ( in the 
case of a separable potential one has exactly A m = and all the a m 's are exactly zero ). The opposite case of A m 
sa Ao corresponds to an accidental situation and is most likely to be found for A m]k and Ao, k belonging to different 
irreducible representations. This should give a smeared second transition below T c which has not been seen up to 
now in high T c superconductors ( except for very recent experiments Jlq] on the penetration depth in YBCO which 
seems to indicate the need of a multiband description; we will come back to this below ). Next we see that in the 
right-hand side of Eq.(8) the logarithm will have a rather small absolute value in most of the range of integration, and 
it will change sign. Moreover in the cases we are interested in, A k also changes sign, and so does A m k in the general 
case. Therefore we have plenty of reasons for destructive interference which will make the right-hand side small in 
general. The conclusion is that taking A k proportional to A , k is a quite good approximation. Moreover replacing 
in Eq.(4) Ao, k by A k should not change much the result when their shape is similar. Hence this should give a very 
good evaluation for the maximum Am of the gap, which is not sensitive anyway to the detailed structure of A k as 
we have seen. 



In order to calculate Am from Eq.(6) we only need to know the weight function N(6) for the reduced gap values 5. 
We introduce the integrated weight x by dx = N(6) dS. In Eq.(6) we can assume by a change of variables that < x 
< 1 . This gives : 

H A M ^dx^xjln^x)) 

nl 1.76T c J /odx^( x ) [ ' 

where 5(x) is a growing function of x with < S < 1 . In the case of two dimensional superconductors, which are 
a very good approximation for all the known high T c superconductors, x is merely the curvilinear abscissa along the 
Fermi line weighted by the local density of states. We can check Eq.(9) in a variety of cases. For a constant gap, S = 
1 and we have naturally the BCS result Am / T c — 1.76 (with 1.76 = tt / e where C is the Euler constant). For the 
A phase of superfluid 3 He, A M / T c = 1.76 e 5 / 6 /2 » 2.029 ( here S(x) = [x(2 - x)] 1 / 2 ). For d-wave, S(x) = sin ( 
ttx/2 ) after change of variable and Am / T c = 1.76 2c" 1 / 2 2.139 . For the simple model introduced by Xu et al. 
||l6|| , with a variable slope at the node, S(x) = x / a for < x < a and S(x) = 1 for a < x < 1 . We find Am / T c 
= 1.76 exp ( a J (9 - 6a) ) which agrees with their result Am / T c = 1.994 for their [i = 4 /(na) — 2 ; for /i = 2.7 
we find Am / T c = 1.904 in agreement with them ; finally for the upper limit a = 1 this gives the highest possible 
value for this model Am / T c = 2.462 . We note that this ratio is not so large, even for this model which does not 
look very physical in this limit ( at least for a single band ). It is then interesting to generalize this model into S(x) 
— x n which gives a very wide gap opening for large n. This leads to Am / T c = 1.76 exp ( n / ( 2n + 1 ) ) , which 
gives 2.631 for n = 2 and saturates at Am / T c = 1.76 exp ( 1 / 2 ) « 2.908 for n — * oo. These few examples show 
that it is quite hard to increase A M / T c even by going to pretty unphysical models. 

In order to explore more fully this question it is convenient to use y = S 2 (x) as a new variable and X(y) = 1 - x as 
a new function. The graph X(y) is trivially related to the graph <5 2 (a:) and it decreases from (y = 0,X=l)to(y 
= , X = 1 ) (even when 5(x) is discontinuous, which occurs for example when the gap is constant ). We can then 
rewrite Eq.(9) into : 
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21n(^) = -l- / > X(y)ln(y) (10) 
[ l-76T c > /o i dyX (y) 

This expression Eq.(10) makes it obvious that, for fixed area J dy X(y) = J dx S 2 (x) , the maximum Am / T c is 
obtained by squeezing as much as possible the weight of X(y) at low y in order to take advantage of the divergence 
of ln(y) for y — > 0. At the same time one sees that this is not very efficient in order to obtain high Am / T c since 
the divergence of ln(y) is weak. This squeezing is optimally reached by taking a constant gap almost everywhere : 
5(x) = 5 m for < x < 1, with 5(0) = and 5(1) = 1, equivalent to X(y) = 1 for < y < 5f n and X(y) = for 5f n 
< x < 1. This conclusion can also be found from Eq.(6) by a convexity argument , as done by Anderson and Morel 
p7[ . The corresponding maximum value of the gap is Am / T c = 1.76 / 5 m . This shows that, by letting 5 m — > 0, 
we can obtain in principle Am / T c as high as we like. However this optimal model is quite unphysical since the gap 
maximum has zero weight, and is therefore irrelevant ( it will not be seen in any experiment ). The real physical gap 
maximum in this model is Am 5 m , not , and we find the BCS value for the gap to T c ratio, which is expected 
since the gap is constant. 

We can consider a slightly more reasonable model by giving a weight 1 - Xo to the gap maximum. In order to obtain 
the optimal Am / T c we take the rest of the gap at a constant value 5 m . Explicitely this leads us to the simple 
model 5(x) = 5 m for < x < xq, and 5(x) = 1 for xq < x < 1. We could try to go continuously from 5 m to 1 , in 
order to obtain a better Am / T c . However it is quite clear that, if we want that the maximum gets a significant 
weight, the improvement will be very small. The above model gives lnr = - xq 5^ ln(<5 m ) / (xo <5„ + 1 - xq ) with 
r = Am / 1.76 T c . We can again in principle obtain r as high as we like by letting 5 m and 1 - xq go to zero. More 
generally, independently of a specific model, it is obvious from Eq.(10) that, in order to obtain a large Am / T c , we 
need qualitatively a gap maximum with a small weight and small gap with a large weight. However the contours of 
constant r of our optimal model are plotted in Fig.l in the ( 5 m , Xq ) plane and they make quite clear quantitatively 
the difficulty which is met when one tries to obtain at the same time a large Am / T c and a sizeable weight 1 - xq for 
the maximum. For fixed r , the maximum possible weight (1 - xo) max is given by (1 - XQ) max = 1/(1 + 2er 2 ln(r)) 
with a corresponding value 1 / (r y/e) for 5 m . While 2 Am / T c = 5 gives a maximum weight 0.2 with 5 m — 0.42 , 
the average experimental value 2 Am / T c = 7 leads to a maximum weight of 0.06 . This is not compatible with the 
experimental data, such as tunneling or Raman scattering which give a gap maximum with a fairly sizeable weight. 
Also the rest of the gap would be at 0.3 the maximum gap value which is rather low. We note that, although the 
above model is already not compatible with experiments, it does not even have nodes in the gap. The following model 
for the gap distribution 5(x) = (x / Xq) 5 m for < x < xq, and 5(x) — 1 for xq < x < 1 , is similar to the preceding 
one, but it is somewhat more realistic since it allows for nodes in the gap. It leads to lnr = x$ 5^ ( 1/3 - ln(5 m )) / 
(xo 5^ + 3(1 - xo))- As expected it gives somewhat worst results for Am / T c . 

Let us summarize the situation. We expect a physically reasonable one-band model to produce a fairly regular gap 
function, similar for example to the standard d-wave order parameter. All the various specific examples of this kind 
that we have considered above gave 2 Am / T c scattered between 4 and 4.5 . Values near 5 correspond already to 
rather unphysical situations. Hence weak coupling is far off the experimental result. Since Am / T c depends on broad 
features of the gap distribution and not on details, as we have discussed, this result is generic not accidental. The 
failure to find higher values than, say, 5 within weak coupling theory is not due to a lack of inspiration in finding the 
proper order parameter. It is a systematic deep limitation of weak coupling theory itself. 

From Eq.(10) the only way to increase Am / T c within weak coupling is to lower the average value of the gap, 
while keeping at the same time a sizeable weight near the gap maximum to obtain agreement with tunneling and 
Raman data. This goes in the direction of a somewhat discontinuous order parameter which does not look like a 
simple one-band model, although one might argue that the spin fluctuation model with strongly peaked interactions 
at wavevector (±7r, ±tt) could produce such a result. We believe that a rather natural realization of such a strange gap 
structure is merely a two-band model (which is included in our study ), with one band corresponding to the maximum 
gap value and the other one to the small value. Nevertheless we have seen that even our optimized model can not 
reproduce at the same time the fairly large Am / T c observed experimentally together with a reasonable weight for 
this gap maximum. We come to the conclusion that simple weak coupling theory is not compatible with experiment. 

What are the ways out ? The most obvious one is to question experiments. As we have seen, this is not an easy 
way since independent experiments are in reasonable agreement. However one may wonder if tunneling or Raman 
experiments do not miss a part of the Fermi surface. In this case the weight of the gap maximum might be less than it 
seems, releasing a part of the theoretical constraint. The most natural situation where this would occur is a two-band 
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model, where the band with the gap maximum would be seen but not so much the other one. On the theoretical 
side one may object that the weak coupling equation Eq.(l) that we have used does not include the possibility of a 
density of states varying strongly perpendicularly to the Fermi surface, as could be produced by nearby Van Hove 
quasi-singularities. However we know that in the isotropic case we have quite generally 2 A m / T c < 4 even for 
such a varying density of states, so that the prospects in this direction are not good. Therefore, if we stop short of 
rejecting BCS theory altogether, the most likely explanation for the high experimental value of Am / T c is that weak 
coupling theory does not apply because strong coupling effects are important. Indeed their existence is supported 
independently by various experiments and they are known to increase in a quite sizeable way this ratio. However 
strong coupling effects with isotropic pairing would have a hard time explaining 2 Am / T c = 7 , since this would 
require jl9[ a quite high coupling constant ( at least 5 ). A possibility is to have strong coupling effects in a multiband 
model. On the other hand self-consistent calculations for simple d-wave pairing within the 2D Hubbard model have 
given results |2(J as high as 2 Am / T c » 10 . Hence it seems that the large experimental Am / T c can be accounted 
for by strong coupling together with anisotropy, whereas it is incompatible with weak coupling theory This comes 
as a strong support in favor of a spin fluctuation mechanism in the debate about the nature of the repulsive pairing 
interaction that we considered in the introduction. 
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FIG. 1. Fig. 1 Contours of constant r = Am / 1-76 T c in the 8 m - xo plane for the model S(x) — 5 m for < x < xq, and 
S(x) = 1 for xq < x < 1. The values of r are 1.1 , 1.2 , 1.3 , 1.5 , 1.7 , 2. and 2.5 as indicated near the curves. 
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